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1. Introduction and motivations 



The strange-quark content of the nucleon is difficult to measure experimentally but is a quan- 
tity of wide interest. In particular, the interaction cross section between some proposed dark mat- 
ter candidates (for instance, neutralinos) and ordinary matter may have a large contribution from 
interactions with sea strange quarks in the nucleon. Specifically, the interest is in the quantity 
(N\J d 3 xss\N) — (0\ J d 3 xss\0): the connected part of the strange quark condensate, integrated 
over the volume of the nucleon. Knowledge of this matrix element is crucial to design experimen- 
tal schemes for dark matter detection and to interpret their results. If the matrix element is known, it 
is possible to determine the constraints on the parameter space of dark matter candidates accessible 
to a given experiment, [g, |||] 

As the quantity in question cannot be measured experimentally or calculated perturbatively, we 
must turn to lattice QCD to compute it. Previous lattice calculations to answer this question have 
been done with quenched simulations^, |[|, with a 2-flavor sea using Wilson or overlap quarks j^, 
[7[ ||], and 2+1 flavor stout quarks^]. Recently this quantity was extracted from baryon mass fits 
to 2+1 flavor simulations JTo|]. The results from these studies often have large uncertainties, and 
some conclude that (N\ss\N) may be significantly larger than its natural size of unity. 

We outline a method for calculating this quantity by evaluating disconnected quark-line dia- 
grams, and present the results of applying this method to the MILC Collaboration's library of gauge 
configurations using 2+1 flavors of Asqtad-improved staggered quarks. An additional very long en- 
semble of 4566 lattices from the UKQCD Collaboration with /3 = 6.75, m\ = 0.06a, m s = 0.30a, 
a ~ 0.125 fm using the same action is included, allowing a measurement of (7Y|jj|/Y) using 25788 
gauge configurations using a full 2+1 flavor sea [{j]]. 

The method outlined here can be used for light quarks as well, and for hadrons other than the 
nucleon. While measurements of their quark condensates are not as immediately in demand, they 
may provide useful information about chiral perturbation theory low-energy constants, and may 
lead to further progress in understanding the structure of the QCD sea within hadrons. 



2. The MILC lattice generation program 

The MILC collaboration is engaged in an extensive project of QCD simulations using a 
Symanzik-improved gauge action and the Asqtad-improved staggered-quark action with a 2+1 
flavor sea. This project consists of a number of runs at nominal lattice spacings of a = .12 fm, 
a = .09 fm, and a = .06 fm, along with other coarser runs not considered here. Details of the 
action, the ensembles of gauge configurations, and the method of extracting nucleon correlators 



can be found in Ref. [11]. An additional very long ensemble of 4566 lattices from UKQCD with 



j3 = 6.75, mi = 0.06a, m s = 0.30a using the same action is included in the analysis. 



3. The method 

In a chiral fermion formulation, the matrix element in question is equal to by the Feynman- 
Hellman theorem. We emphasize that the derivative should be taken holding all other parameters 
in the action fixed. A sizeable value for does not imply a large dependence of the physical 



2 



2+1 flavor nucleon strangeness 



Walter Freeman 



nucleon mass on m s . Changing the value of m s changes the value of all physical quantities by a 
similar amount, and this change is interpreted as an overall rescaling of the lattice rather than a shift 
in the physical values. Likewise, a large result does not imply the presence of many ss loops in the 
nucleon. Rather, it comes from the suppression of the vacuum ss condensate near the nucleon. 

One needs large-scale simulations with a true 2+1 flavor sea to determine (iV|jj|iV) with rea- 
sonable accuracy. Previous methods used to determine this quantity require particular choices for 
the lattice parameters, namely ensembles with different m s but the same j3. On the other hand, 
in a general program of lattice QCD simulations such as MILC's, m s is typically held fixed at its 
physical value while /3 is changed with the m 9 's to keep the physical lattice spacing fixed. Thus, 
computing (iV|,M|iV) using the standard method requires special-purpose simulations which must 
be limited in scope due to economics. We have developed a method to determine from any 
single lattice ensemble, allowing its use on the MILC lattice configurations. The many ensembles 
are simply used to conduct an extrapolation to the physical point and to improve statistics. 

The nucleon mass Mn is obtained by a fit to the nucleon correlator C(t) and as such can be 
thought of as a complicated function of the correlator at different times: Mn = f(C(h ) , Cfa) , Cfo) . . .) 
The crucial idea is that one can use the chain rule for differentiation to rewrite the derivative: 



dM N _ dM N dC(h) | dM N dC(t 2 ) | dM N dC(t 3 ) 
dnig dC(t\) dm s dCfa) dm s dC(t$) dm s ... 

The partial derivatives can be evaluated most simply by applying a small perturbation 
to the nucleon correlator and examining the change in the fit result. The other partial derivative 
can be evaluated by an application of the Feynman-Hellman theorem in reverse to relate it to 

(p m (t i )ss)-(P(t i ))(ss). 

Here we take advantage of the fact that, whenever the MILC code generates or reads a lattice 
for analysis, it prints a stochastic estimator for J d A xss. The number of estimators used in this 
work, 6-16 per lattice, is sufficient for their fluctuation to contribute significantly to the statistical 
error. These values can be used to compute (C(t) J d 4 xss) — (C(t)) (/ d A xss) by simple evaluation 
with no additional use of computer time. By doing this for each t used in the fit to determine Mn, 
it is possible to evaluate the chain rule sum and determine 

This double use of the Feynman-Hellman theorem to access the nuclear strangeness might 
seem redundant: as we have the values of (ss), why not simply evaluate the matrix element (N\ss \N) 
directly? The problem is that the lattice operator used to create and annihilate the nucleon overlaps 
with many other three-quark states; the nucleon is simply the lowest-lying three-quark state. Fur- 
thermore, the normalization of the nucleon state created this way is unknown. Using this operator 
along with the (ss) data to evaluate (A^|5J|A^) directly would give the sea quark content of some 
superposition of states of unknown normalization, not of the nucleon itself. The fitting procedure 
used to extract the nucleon mass provides a way of extracting information about the nucleon alone, 
and the double use of the Feynman-Hellman theorem provides a way to relate it to (N\ss\N). 

Statistical errors on the result for were calculated using the jackknife method with blocks 
of size 10; this is large enough to take into account autocorrelations in the lattice data. Use of larger 
jackknife blocks produces an insignificant change in the size of the errors. 
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4. Choice of fit range 



Choosing a lower minimum distance in the nucleon mass fits will result in lower statistical 
error but may introduce systematic biases from pollution by excited states. The minimum fit dis- 
tance required to avoid significant systematic error in (iV|jj|iV) can be smaller than the one used for 
precision measurements of M#. When extracting M#, any pollution of the correlator by the excited 
states in the fit range will cause an incorrectly high value for the mass. However, the lowest-lying 
excited state is the delta. As we expect the effect of the delta on the strange quark condensate to be 
broadly similar to that of the nucleon, a small amount of delta pollution will not create substantial 
systematic error. 

The minimum distances chosen should be consistent in physical units between lattice spacings. 
The minimum distance chosen should also be in a region where the result at successive minimum 
distances does not differ greatly. Both looking at the result ensemble-by-ensemble and looking at 
fits to all ensembles in the same nominal lattice spacing suggest the use of t m i„ = 0.6 fm. This 
choice is also suggested by the signal-to-noise ratio in T^p; see Figure [j]. In lattice units, this 
gives t m i„ = 5a (a = 0.13fm), t m i„ = la{a = 0.09fm), and t m {„ = 10a (a = 0.06fm). By examining 
the dependence of on f min , we conservatively estimate the systematic error due to excited state 
pollution as 10%. However, the nucleon mass itself can be computed with much lower statistical 
error at higher ? m; „; mass fits at these higher minimum distances differ from those at t m i„ = 0.6 fm 
by only 1% to 5%. Thus 10% is potentially an overly-pessimistic estimate. 

The result is quite insensitive to the maximum distance used for the fits; this is expected, since 
the signal-to-noise ratio of the correlator is very poor there. 

The values of obtained on each ensemble using these minimum distances are tabulated 

am s ° 

here and shown in the first panel of Figure ^[ In the table, both the bare value (in the lattice 
regularization), and the value converted to the ms(2 GeV) regularization and shifted to the correct 
strange quark mass (see below) are given. 

T 



~6xr 

6.79 
6.76 
6.76 
6.76 
6.75 
~7M~ 
7.085 
7.09 
7.10 
7.11 
7.18 
7.10 
~7M~ 
7.465 
7.47 
7.475 
7.48 
7.46 



ami 
0.30 
0.20 
0.10 
0.07 
0.05 
0.06 
0.0031 
0.00465 
0.0062 
0.0093 
0.0124 
0.0310 
0.0062 
0.0018 
0.0025 
0.0036 
0.0054 
0.0072 
0.0036 



am s 
0.50 
0.50 
0.50 
0.50 
0.50 
0.30 
0.031 
0.031 
0.031 
0.031 
0.031 
0.031 
0.0186 
0.018 
0.018 
0.018 
0.018 
0.018 
0.0108 



a (fm) 
0.117 
0.118 
0.119 
0.118 
0.117 
0.117 
0.084 
0.084 
0.084 
0.084 
0.084 
0.081 
0.082 
0.059 
0.059 
0.058 
0.059 
0.059 
0.058 



Nlats 

~55T~ 
2067 
2278 
2098 
2033 
4566 
~T0T3~ 
599 
1943 
1137 
1993 
496 
948 
STT 
798 
658 
616 
620 
550 



(NssN) (bare) 
0.676(190) 
0.702(98) 
0.779(137) 
0.867(214) 
0.753(299) 
0.884(171) 
1.232(339) 
0.500(369) 
0.705(158) 
1.093(183) 
0.936(109) 
0.530(161) 
0.776(218) 



0.556(366) 
0.848(494) 
0.561(268) 
0.813(443) 
1.352(249) 
0.643(438) 



(MssN) (adj) 
0.620 147) 
0.639(76) 
0.696(106) 
0.766(166) 
0.679(230) 
0.645(132 
0.955(249 
0.417(271 
0.568(116 
0.853(134) 
0.739(80) 
0.449(118) 
0.500(159) 
0.375(255 



0.579(345 
0.379(187) 
0.554(309) 
0.929(160) 
0.328(305) 



5. Analysis 

5.1 Renormalization 

Since in the continuum extrapolation it is necessary to compare quantities measured at different 



4 



2+1 flavor nucleoli strangeness 



Walter Freeman 




Figure 1: The nucleon correlator and the derivative of this correlator with respect to m s for the ensemble with 
ami = 0.0093 and am s = 0.031 (first panel). For the derivative, the squares are points where the derivative 
is negative, and crosses are points where it is positive. The vertical lines show the range used in fitting the 
correlator. The second panel shows for three ensembles with a w 0.9fm as a function of the minimum 
distance used in the fitting, and the third panel shows the fitted nucleon mass itself versus ?„„•„. The error 
bars labelled "10%" in the second and third panels show the size of the ten percent systematic error estimate 
from excited state contamination. 




Figure 2: -j^- evaluated on the MILC ensembles. The first panel shows the "bare" results; the second 
panel shows the results converted to the ms(2 GeV) renormalization and adjusted to the physical strange 
quark mass. The dotted line shows the chiral and continuum fit; the black point at m; = .005 shows this fit 
evaluated at the physical point. 



lattice spacings and thus different regularization schemes, it is necessary to convert all of the data 
to a common renormalization scheme since ^ is a renormalization dependent quantity. Since the 
final result will be presented in the MS renormalization scheme at a scale of 2 GeV so it will be 
most useful to the broader physics community, it is preferable to begin the analysis by converting 



all values to this renormalization scale. The Z-f actors for this conversion are known [|12|]. 

5.2 Strange quark mass dependence 

Since the value of m s in lattice units is not known prior to a lattice run, all of the ensembles 
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were run at strange quark masses different than the physical one; on some ensembles this error 



is large (20%). It is thus necessary to determine the dependence of -f^- on m s (that is, to 



perform an extrapolation to the physical m s . 

As MILC has run one or two ensembles at each nominal lattice spacing with m s set to 60% of 
the guess at the physical value, it should be possible to determine by examining on both 



these ensembles and those with the heavier strange quarks. However, since these "light strange" 
ensembles are short, better results for this second derivative can be obtained by a different method. 

On many ensembles, the light quark is quite heavy; in some cases, it is actually closer to the 
physical value of m s than the lattice heavy quark. Thus the behavior of the "light" quark condensate 
on these ensembles is similar to that of a strange quark condensate with a lower strange quark mass. 

These ensembles with a heavier light quark allow for the measurement of the behavior of 
a "heavy quark condensate" at two different quark masses, and allow us to estimate the second 
derivative ^fy- = - l(m s - m,). Since ^ and are measured on the same lattices, 

amj K om s ami ' 1 v J 1 ' am s ami ' 

their values are correlated, so the error on — (determined by jackknife) is reduced. 

By examining five ensembles with a large number of configurations and relatively heavy light 
quarks, it is possible to determine d ^ = —2.2 by a weighted average of its value on each, esti- 
mated as above. This value is then used to extrapolate each data point to m SiP hy S - The second panel 
of Figure |2| shows the data converted to MS(2GeV) and adjusted to the correct m s . 

5.3 Light quark mass dependence and continuum extrapolation 

Similarly, the value of may depend on mi, and we are most interested in evaluating it at 
the physical light quark mass. In this case the data set contains results for ^ L at many different 
values of mi, so it is possible to determine the dependence on mi with a simple fit. Examination 
of the #PT form for reveals a constant plus linear fit in m; is good enough; no #PT terms at 



higher-order are relevant at the level of statistical accuracy provided by the present data [13]. 

It is also necessary to extrapolate to the continuum. In the Asqtad fermion formulation, the 
leading-order errors in the action are proportional to a 2 . Thus, the leading-order effect on 
will likewise be proportional to a 2 ; this effect can be determined by adding such a term to the fit 
form. However, such a term will be poorly constrained, since the better statistics in our dataset are 
from the a = .12 fm ensembles, and since the effect from lattice spacing is small. Other hadronic 
quantities calculated using the MILC Asqtad data show roughly a 10% effect between the coarse 
lattices and the continuum, so we use a Bayesian prior with a central value of and a width 
corresponding to a 10% effect to constrain the lattice spacing dependent term in the fit. 



As both the light quark mass dependence and the lattice spacing dependence will be computed 



in the same fit, the proper fit form is = A + 6m/ + Ca 2 , with a Bayesian constraint on C. 



6. Result and error budget 

Evaluating the fit above at the physical value of mi and in the continuum, we find = 
0.69 ±0.07 statistical- We estimate the systematic error due to excited states present in the nucleon 
correlator, as discussed above, at 10%. The extrapolation to the physical light quark mass involves 
higher order terms in %PT which were not considered here. To estimate the size of this effect, we 
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note that if the nucleon mass itself is fit to a constant-plus-linear form over the range considered 
here, the result is seven percent off from the result obtained when two more orders in m n are added 
to the fit. We thus estimate the effect of higher order terms in #PT as 7%. In one case where 
a spatial volume larger than the one used here has been run, the nucleon mass computed on the 
larger volume differs by 1%. Since the effect on the strange quark condensate is potentially more 
sensitive to finite volume effects, we estimate the systematic error due to finite volume effects as 
3%. Finally, Ref. [|l2|] quotes an error in Z m as 4%. If these errors are combined in quadrature we 
thus estimate the total systematic error as 0.09. 

The renormalization-invariant quantity m s ^f- is also commonly quoted. Using a similar fit, 
we calculate m s ^^- = 59(6) (8) MeV. This quantity does not include uncertainty this 
cancels, but includes a lattice systematics error of nearly the same size, coming from the 2% un- 
certainty in the lattice scale and the uncertainty in the lattice strange quark mass. 
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